Nonequilibrium Multi-Band Spin Quantum Transport Equations: 
Spin, Pseudo-Spin, and Total Charge Coupling 

Felix A. BuotQ 

Center for Computational Materials Science 
Ceorge Mason University, Fairfax, VA 22030, U.S. A 

and 

C&LB Research Institute, Carmen, Cebu, Philippines 

Abstract 

Using the superfield nonequilibrium Green's function technique, we derive the spatio-temporal 
spin magnetization quantum transport equations (SMQTEs) for a two-band model of semicon- 
ductors. The relevant variables are the real (Pauli-Dirac) spin, pseudo-spin, and the total charge. 
The results show that the multi-band real SMQTEs are coupled to the pseudo-spin magnetization 
transport equations by virtue of the presence of two additional discrete quantum labels besides 
the up and down real-spin indices, namely, the conduction and valence band quantum labels. The 
SMQTEs essentially consist of three group of terms describing the rate of change, namely, (1) a 
group of terms similar to the equation for particle quantum transport, i.e., with spin-independent 
transport parameters, (2) a group of terms describing various torques influencing the spin orienta- 
tion and directional flow of spin magnetization correlations or phase-space magnetization density, 
and (3) a group of terms expressing the coupling of the real spin magnetization with the pseudo- 
spin magnetization. Self-consistently, the pseudo-spin magnetization equations incorporate the 
pseudo-spin/real spin coupling, as well as the pseudo-spin coupling to the total charge. 
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I. INTRODUCTION 



The need to analyze the ultrafast-switching-speed and power-dissipation (speed-power 
product) performance of nanoelectronic devices has ushered the extension of the classical 
Boltzmann transport equation to a fully time-dependent and highly-nonlinear nonequilib- 
rium quantum distribution function (QDF) transport equations for charge carriers. This has 
been achieved through the use of non-equilibrium Green's function, obtained either by the 
time-contour quantum field formulation of Schwinger^, Keldysb^, and Kadanoff and Baynt^, 
or by the real-time quantum superfield formulation of Buot, coupled with his lattice Weyl 
(LW) transformation technique.-"- This extension to QDF transport equation has proved to 
be highly crucial in discovering autonomous THz current oscillations in resonant tunneling 
devices through numerical simulations, and in resolving controversial issues concerning the 
highly-nonlinear and bistable current-voltage characteristics found in the experiments^*^. In- 
deed, in the phase space QDF kinetic approach has so far been the most successful technique 
in the time-dependent analyses of open and active nanosystem and nanodevices, as evidenced 
by the numerical work of Jensen and Buot on resonant tunneling heterostructures^. 

To the authors' knowledge, the extension of the classical spin density equation or the 
Bloch equation for spin transport^, which is the analogue of the classical Boltzmann equa- 
tion for charged-particle transport, to fully space-time dependent and highly-nonlinear 
nonequilibrium QDF transport equations for the magnetization, has only been recently 
reported by the author for a single energy band in a co-authored paper— . With the ex- 
ploding surge of interest on spintronics and nanomagneticsi^, there is an urgent need for 
this fully quantum transport extension of the classical Bloch equation to guide the time- 
dependent numerical simulation of the speed/power switching performance and reliability 
analyses of realistic functional spin nanostructures and transistors. In this paper, we report 
the multiband SMQTEs. 

An immense activity with time-dependent spin magnetization equation, driven by techno- 
logical applications, occurs in the field of micromagnetics dealing mostly with micro-domain 
wall dynamics. The so-called classical Landau-Lifsitz-Gilbert equation, with its various 
modifications and phenomenological damping terms, has become the major player in the 
analyses of the dynamics of the magnetization field.—"— 

In magnetoelectronics and spintronics, the mutual dependence of electron transport prop- 
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erties and magnetic properties have resulted in various phenomena referred to as giant 
magnetoresistancei^, spin modulator— based on spin precession with or without magnetic 
field (based on spin-orbit coupling), and spin transfer torquei^*^ which have become spring- 
boards of several novel spin-based device concepts. 

In analyzing the current-induced spin orientation of electrons in semiconductors, with 
spin-orbit coupling, Dyakonov and PereU^ used the following transport equation for the 
spin spatial- density vector, S, as 



dt 



-v-g — 

Ts 



+ nx s, (1) 



+ -fM X B, where Mr = <—,^,— } , 2 

I T2 T2 Ti J 



where Q is the spin flux density, 

n 

/is is the Bohr magneton, g is the g-factor for electrons, and B is the effective magnetic 
field. Earlier, in the 1950 's Torrey^U employed the classical Bloch equations for spin systems 
in the diffusive regime as, 

dM 

where D is the diffusion coefficient, ti is the spin- magnetization relaxation characteristic 
time, and T2 is the spin-dephasing characteristic time for the M^ and My components. Most 
of the works that follows on spin transport also make use of classical spin transport equations. 
These classical techniques are no longer valid for analyses of the ultrafast switching speed 
and power-dissipation performance of the emerging spintronic devices. 

What these classical treatments have earlier shown is that spin transport consist of 
terms similar to particle transport, i.e., with spin-independent transport parameters, such 
as the diffusion coefficient, D, and relaxation times in Eq. ([2]), and terms describing the 
torques in the system.— However, the spin-charge coupling for a single band transport^, 
and spin/pseudo-spin/charge interaction^^ in multiband transport are lacking in the equa- 
tions or at best not fully treated self-consistently within these classical and semiclassical 
treatments.— 

The purpose of this paper is to extend the above classical magnetization transport 
equations to a fully time-dependent and highly-nonlinear nonequilibrium spin magnetiza- 
tion QDF transport equations in phase space, based on nonequilibrium Green's function 



technique,-"-'^ on the same level as in the extension of the classical Boltzmann kinetic equa- 
tion to full QDF transport equation for electrons. In the single-band magnetization QDF 
transport equations, the spin-charge coupling is explicitly incorporated. In our present 
multi-band generalization, we found that besides the highly-coupled transport equations for 
the Pauli-Dirac or real-spin density, nonlinear coupling to the pseudo-spin density is also 
present in all of the spin magnetization quantum transport equations for coupled electrons 
and holes. Conversely, the pseudo-spin magnetization quantum transport equations incorpo- 
rate the pseudo-spin/real-spin coupling, as well as the pseudo-spin coupling with the total 
charge in the system. The SMQTEs are summarily cast in more physically transparent 
forms in Sec. IVIIII 



II. NONEQUILIBRIUM QUANTUM SUPERFIELD THEORY 

Our starting point is the general quantum transport expressions for fermions as obtained 
from the real-time quantum superfield theoretical formulation of Buot^'^: 

+ [Single - glK^l] ■ (3) 

The last two brackets account for the Cooper pairings between Fermions of the same specie. 
These do not concern us in this paper (their corresponding transport equationa^"^ are im- 
portant in nonequilibrium superconductivity). In what follows we will drop these last two 
brackets of the RHS of Eq. 
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III. MULTI-BAND QUANTUM TRANSPORT EQUATIONS 



In the absence of pairing between fermions of the same-specie Eq. ([3]) becomes, by 
exphcitly writing the quantum arguments, as 

= [vJ^,6,^G% (22) - G% (12) <5^^<522] 

+ [s;^ (12) G% (22) - G% (12) s:;, (22)] 

+ (12) G:;^ (22) - Gl^ (12) (22)] , (4) 

where the Greek subscript indices correspond to discrete band indices, and the numeral 
indices correspond to the two-point space-time arguments. In what follows we will treat 
the two-band model of a semiconductor and replace by v and c, the Greek indices for the 
valence and conduction band quantum labels, respectively. 



A. Electron-Hole Picture 



The transport equation for the conduction-band electrons in the electron-hole or defect 
representation reduces to 



^dti dt2 



Glt^< (12) 

= w (10 G'< (^2) - G< (10 vi m\ 
+ (iOG<(e2)-G<(io Kcim 

+ [E< (10 Gl, m - (10 S< (^2)] 

+ (10 (^2) - ^-t< (10 K-eit m 

+ AM,t< (10 ^^Ltr (^2) - g^^-^^ (10 K-elt m 



(5) 



where the off-diagonal self-energies, S^^, and Green's functions, G^/j are denoted by 
and f/^"^, respectively, in analogy to the A-function and anomalous Green's function, J-", 
i.e., g^"^ ^ J-" in the theory of superconductivity.—"— In Eq. (E]), there is no Cooper 
pairing between electrons; moreover, the electron and electron-hole pictures coincides for 
the conduction band. 
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To obtain the equation for the hole density from the general superfield formulation for 
electrons, we make use of the electron-hole conversion table given in Ref.-i^, where the 
relevant portion is reproduced below, 



Table 1. Mapping from Electron to Electron-Hole Picture 



electron picture 


(e — field) 


{e — h field) 


e — h picture 


-*^G'< (12) 


(^i (2)^.(1)) 


(0.(2)0Ul)> 


^fiGt>^ (12) 


-znG< (12) 


(V'J (2)^.(1)) 


(^cU2)<^Ul)> 


-ihgl-;:;^< (12) 


-zhG< (12) 


(^i(2)^c(l)> 


(0, {2)ij,{l)) 


(12) 


-znG< (12) 


(V'K2)^e(l)> 


(^l(2)^e(l)> 


-ihGl-^^< (12) 


^^G>J(12) 


(V^.(2)V'i(l)> 


(0i(2)0.(l)> 


-z^G„^;< (12) 


z^G> (12) 


(^.(i)V'i (2)> 


(0Ul)0. (2)> 


-ihG'if^ (12) 



We also have the following Tables, which can also be similarly applied to the self-energies. 



Table 2. 



Table 3. 



electron picture 


e — h picture 


Gl, (12) 


-Gl-^-'^^ (12) 


Gl, (12) 


9ltvc (12) 


(12) 


i^M,c;"(i2) 


Gl. (12) 


GL-''M12) 



electron picture 


e — h picture 


Gl, (12) 


-g:-^'^^ (12) 


G."c (12) 


^7:etcMi2) 


Gt. (12) 


^/h^;" (12) 


G'?c(12) 


^^-^^'"(12) 



From the first table, we see that the hole density, —iHGl'^'^''^ (12) correspond to ifiG^J (12) 
in the electron energy-band representation. Thus, to obtain the equation for the hole den- 
sity given by —iHGl~^'^ (12) , one may look for the equation for ifiG^J (12) in the general 
nonequilibrium formulation of electrons. The equation for iHG^J (12) can be obtained from 
the equation for G^^,by taking the complex conjugate of the above equation and making use 
of the relation 



G>1-(12) = -G>'^'(12) 



(6) 
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We have the transport equation for (12) given by 

= [vl (10 G'> (^2) - G> (10 v.. m] 
+ [S^, (10 G> (^2) - G> (10 (^2)] 

+ [s;,(iOG> (e2)-G> (ios^.(e2)] 

+ [s> (10 Gi m - g:^ (10 s> (^2)] 

+ [s> (10 g:, m - g:, (10 s> (^2)] (7) 

Taking the complex conjugate of Eq. ([7]), we obtain 

= [-G> (20 vl (O) + ^J. (^2) G>J (10] 
+ [-G'> (20 + (20 G>, m 

+ [-G> (20 (ei) + (20 Gl m 
+ i-Gl, (20 S> (O) + S> (20 Gl m 

+ (20 S> (ei) + S> (20 Gl (O)] . (8) 

Applying the relation 

-G>^ (12) = -G> (21) = G-^-< (12) , 

and going entirely to the defect representation for the rest of the terms, we obtain for the 
time evolution of the hole density without Cooper pairing between holes. 

= - (10 g:-'^'< m + Gi-^'< (le) vi m] 

- [-^-^^^ (10 Gi-^< m + Gi-^< (10 K-.^'" m] 

- [-T^i-^'^ (10 Gi-^^< m + Gi-^< (10 k;;^'" m] 

- [-Kt'< (10 Gi-^^- m + Gt^^' (10 m] 

- [-S^-'^'< (10 Gl-^^'^ m + Gl-^^^ (10 S:-'^'< (^2)] . (9) 
Replacing the off-diagonal self-energies, S^^, and Green's functions, G^p by A"^"^ and (/'^~'^, 
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respectively, as was done in Eq. ([5]), allow us to rewrite the equation as 

- [E^;'^'^ (10 Gi-'^^< m - g:-'^'< (10 K^''^ m] 



m 9ii%^ m - siif (10 Ae 



vv 

hh,cv 



-h,a 



[s:;'^'< (10 Gi-'''^ m - g:-'^'^ (10 s:;'^'< m] 



,e—h,< 
ee,vc 



(10 gii!^: m - 9::X (lo k;;'"^ m 



. e~h,< 



(10) 



Compared with Eq. ^ for the conduction band, one can readily see that for fiat bands, or 
atomic limit, the residual terms in these two equations are equal. 



dti dt2 



Gl 



Gl 



'h,< 



-h,< 



(12) 



(12) 



\e—h,r e—h,< 
ee'uc yhh,cv 



e-h,<A e-h,a 
iiee,vc hh,cv 



e-h,r e-h,< 
hh,cv yee,vc 



e-h,< 
9hh 



/\e—h,a 

cv ee,vc 



A 



e—h,r e—h,< 
hh,cv yee,vc 



e—h,< \e—h,a 
yhh,cv ee,vc 



(11) 



Since A^,-V = A^,-^^^ and Ae^t^ 



^hhcv- This equality simply states that the rate of 



change of the holes in the valence band is the same as the rate of change of the electrons in 
the conduction band, or the rate of electron creation in conduction band equals the rate of 
hole creation in the valence band, and vice versa. 

If we use the electron picture to write the transport equation for the valence band, the 
resulting equation would be similar to Eq. ffTOj) without the superscript ^"'^ and with the 
right hand side of that equation multiplied by —1. Then one can also easily see that for 
fiat-band case, the rate of change of electrons in the conduction band and that of the valence 
band would have opposite sign, i.e., the rate of creation of electrons in the conduction band 
occurs at the expense of equal rate of loss of electrons in the valence band. Later in this 
paper, we will use the electron picture to investigate the pseudo-spin transport equations. 
There, the total charge represents the total net charge to be used in the Poisson equation. 
Of course the total net charge at each point of the phase-space is the charge represented by 
the total number of electrons in the conduction and valence band minus the total positive 
background charge of the crystal lattice. 

Note that the single particle Hamiltonian v^^ is not really a two-point function, unlilce 
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the correlation functions. The above equation reduces to 



+ [E-'^'^ (10 Gi-^^< m - Gi-'^'< (10 Kv"'"" m] 
+ [Ei-''< (10 Gi-^'^ m - Gi-^^^ (10 



+ 
+ 



(20 Kr!:^ m 



,vc hh.cv 



(20 gJui: (a) 



A:-.'c>(20^^i;"(ei)]- 



(12) 



The equations for the interband or polarization terms, gl^^'^i ^tc. can be obtained 

from the equations for (12), where a 13. In the electron-hole or defect representation, 
the equation for g^^^^''^^ is determined from the general superfield formalism given by the 
equation for G^, (12). Upon transforming to the electron-hole picture, the change of polar- 
ization due to the destruction of electron-hole pairs brought about by Auger recombination, 
optical de-excitation, and other recombination processes is given in the defect representation 
as. 



ih 

+ 
+ 
+ 
+ 



e-h,< 
9 hh.cv 



(12) 



m 9hK^ (10 - m vi (10 
Kc''- (10 fhiv^ m - Gi;''< (10 Kh%^ m 
^17''^ (10 gVht^ m - Gi7''^ (10 Ki.^ m 



-KhX (10 Gi-''>^ m + (10 Kv'^"-^ m 



-Kh^ (10 Gl-''^"^ m + (10 E-'^'>^ (^2) 



(13) 



Similarly, upon going over to the defect representation of the reverse process, i.e., the 
change of polarization due to the creation of electron-hole pairs by impact ionization, optical 
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excitation, Zener tunneling, and other excitation processes is therefore given by 

+ i-K^''^^ m 9i:if m + Gi-'^'>^ (10 K-eif m] 
+ [-K-j'^>^ (10 gt^i^ m + Gi-'^''^^ (10 Ktvf m] 
+ [A-Y (10 g:-^^ m - ^-^^< (10 Kc^^" m] 

+ [a:-::;< (10 Gi:;:^^'^ m - (10 s^,-'^'< (?2)] . (14) 

Note the presence of the 'transposed' terms, namely, G^-^'^^ (^2), S^;'^'''^ (^2), Gl^^'""^ (^2), 
and E-'^'>^ (^2) in Eq. m and S^'^''^^ (10, Gl'^^'^^ (10, G-'^'>^ (10, and S-'^'>^ (10 
in Eq. ( IT^ of the interband 'pairing' correlation functions. 

In the calculation of the pseudo-spin transport equations the transposed quantities in 
Eqs. ( IT3l) and ( 1T4|) do not enter since we use the electron picture to provide the two states 
for the electrons, namely the two quantum labels c and v in order to derive the pseudo-spin 
quantum transport equations. Since there is a one-to-one mapping between the electron 
picture and electron-hole picture except for the valence band as shown in the tables above, 
all one does to go from electron-hole picture to the electron picture is to change the sign of 
the right-hand side of Eq. ( |T2|) and substitute for the transposed quantities their equivalent 
expressions in the electron picture. 

IV. EQUATIONS FOR BLOCK ELECTRONS WITH SPIN 

In the presence of Pauli-Dirac spin degree of freedom, Eqs. ([5]), ( fT2|) . ( 1T3|) . and ( IT^ 
become matrix equations. These are given in the Appendix. 

A. Spin Canonical 2x2 Matrix Forms 

In spintronics, we are interested in the time-dependent evolution of the multi-band polar- 
ization and magnetization densities, Sap^z = ifi {G^f^ ^-f " G'^^ , as these are transported 
across the device. This leads us to transform all 2 x 2 matrices into their spin canonical 
forms, defined here as expansion in terms of the Pauli spin matrices and identity. 
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First let us partition the total nonequilibrium Green's functions (TNEGF) into 2x2 
submatrix components, by virtue of the spin indices, as 




(CC) (CV) 

(yc) (vv) 



(15) 



where (CC), (CV), (VC), and (VV) are 2x2 submatrices. We can form four coupled 
spin magnetization quantum transport equations for each submatrix, giving us 16 coupled 
magnetization quantum transport equations. 

Thus, we first transform the 2x2 matrix variables in Eq. f ll5p into their spin-canonical 
forms using the Pauli-matrices, 



where we dropped the superscript ^ in the above spin-canonical form, this is to be understood 
in what follows unless otherwise specifically specified. The spin canonical expansion in terms 
of the Pauli matrices basically separates the spin-independent terms from the spin- dependent 
terms. 

In the case of the conduction, Eq. (fT6l) . and valence bands, Eq. (fT9|) . the coefficients of 
the identity matrix, namely, Scc,o and Svv,o represent the charge density of electrons in the 
conduction band and charge density of holes in the valence band, respectively. Similarly, 
Scv,o and Svc,o represent the probability density for the annihilation and creation of electron- 
hole pairs, respectively. The coefficients of the vector, a, represent the corresponding spin 
magnetization vectors. Thus the canonical form basically performs the essential function of 
separating particle charge and spin. 

We will later show that the coefficients of the identity matrix given above can also be 
expressed in terms of the pseudo-spin vector and the total charge in our multi-band system. 




(19) 



(17) 



(18) 



(16) 



11 



B. Intraband Nonequilibrium Spin Correlation Functions 



To establish the notations used in this paper, we give the spin canonical forms for the 
single particle Hamiltonian, various correlation functions, and self-energies. 

For the intraband canonical terms, we have for the intraband spin-subtraces scalar cor- 
relation functions, 

•^c^^^^ = {g:-'^^'''^< + G:-if'''^<) , (20) 

Qr,a,< _ ( /~ie-h,r,a,< . ^e-h,r,a,<\ /q-. \ 

The intraband spin-correlation vector components for electrons are, 

qr,a,< _ ( ^e-h,r,a,< , ^e-h,r,a,<\ 

■ Qr.a,< I f^e—h,r,a,< ^e—h,r,a,< 



and for holes. 



qr,a,< _ ( y^e—h,r,a,< _ ^e-h,r,a,<\ 

'^cc,z — y^cc,n ^cc,U ) ' 

cr,a,< _ f ^e-h,r,a,< , ^e-h,r,a,<\ 

■ qr,a,< / f^e—h,r,a,< ^e—h,r,a,< 

qr,a,< / ^e-h,r,a,< _ y^e-h,r,a,<\ 



C. Interband Nonequilibrium Spin Correlation Functions 

For the interband canonical terms, we have interband spin-subtraces scalar correlation 
functions, 

qr,a,< _ f e-h,r,a,< , e-h,r,a,<\ (r,r,\ 
'^cv,o — \ yhh,cv,r[ ~'~ yhh,cv,U J ' ^ ^ 

qr,a,< _ f e-h,r,a,< , e-h,r,a,<\ /r)q\ 
'^vc,o — \yee,vc,U yee,i;c,tt / ' ^ ^ 

and for the interband spin- vector correlation components, 

qr,a,< _ f e-h,r,a,< , e-h,r,a,<\ 
— \yhh,cv,X't ~'~ yhh,cv,U J ' 

•cr,a,< _ f e-h-r,a,< _ e~h,r,a,<\ 
'"^cv,y — yyhh,cv,i^ yhh,cv,U J ' 

qr,a,< _ f e-h,r,a,< _ e-h,r,a,<\ 

'^cv,z — yyhh,cv,n yhh,cv,u ) ■ 
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The corresponding interband reverse process 'conjugate' correlations are, 



cr,a,< 

vc,x 


I e-h,r,a,< 

I "66,1)0, t4- 


+ 


■qr,a,< _ 
vc,y 


f e-h,r,a,< 


-9, 


qr,a,< _ 
vc,z 


f e-h,r,a,< 
\yee,vc,].\. 


-9, 



tt J ■ 



D. Single-Particle Hamiltonian Spin Canonical Forms 

We have for the single-particle Hamiltonian for the electrons expressed in canonical forms, 
assuming the presence of external magnetic field and/or spin-orbit coupling. 



where the upper bar in indicates the subtrace of the 2x2 spin matrix for the conduction 
band. Similarly, for the holes we have, 

Vv,aa' = ^ (^H^I + By-a^ . 
These lead to the matrix expression for electrons. 



He ~\~ Bc^z ^c,x ^^c,y 
Vc,(T(t' ~ 2 \ - ' ' 



and for holes. 



1 / Hy -\- }3y^z ^V,X iBy^y 

Vv,aa' = 2 I 

The components of Be and B^ are real valued, since the single-particle Hamiltonian, v, 
is Hermitian. These are the effective magnetic field components, multiplied by HbQi which 
accounts for the external magnetic fields as well as the effects of spin-orbit coupling. 

E. Canonical Forms for the Spin-Dependent Electrons Self-Energy 

Similarly, we have for the spin-canonical forms for the electron self-energies, 

1 / _(_ '^r _ j'^r 
Y^r j cc ~ "cc,z ^cc,x cc,y 

2 I '^r _|_ j'^r -yr '^r 

\ "cc,x ' cc,y cc "cc,z 



cc.aa' 
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ya _ 

cc,aa' 2 



1 / +s 



cc ' — cc,z — cc,x " — cc,y 
■"^cCjO: "I" cc,y cc ^cc,z 



Yi< cc ^ "cc,z "cc,a; cc,y 

cc.aa' o I , ^ - ^ 

\ "cCjO: ~ cc,y cc ^cc,z 



F. Spin-Canonical Forms for the Holes Self-Energy 



vv,aa' 2 



vv ' — vv,z — vv,x — vv,y 
vv,x "I" vv,y vv " till, 2 



1 / y« -I- 

Y^a I I'll ^ " ,2 " vv,x 1)11,3/ 

vv.aa' o I - 

\ ^ vv,x ' vv,y vv " vv,z 



1 / y< -I- " < "< _ 7 "< 

■r-i< ( VV ' " vv,z " vv,x vv,y 

vv,aa' o I ^ ^ - ^ 

^ DD,a; ' 'U'y,2/ uf ' — ' vv,z 



G. Canonical Forms of Electron-Hole Pairing Self-Energy 



^hh,cv,aa' ~ ^ \;^hh,cv ^ ^hh,cv ' ^ j i 

^hh,cv,aa' ~ ^ {^hh,cv + ^hh,cv ' ^) ' 

\e—h,a ^ / X a i ?a _ 

ee,vc,(7(7' 2 \^ ee,i;c ~r "ee,i;c ' ^ J ' 

A e-/j,< ^ 1 /^A< I ?< . ^^ 
ee,?;c,CTt 2 \ ee,i)c ' "ee,vc ^ J ' 



\e—h,r j hh,cv ' hh,cv,z hh,cv,x hh,cv,y 

A?" _|_ j^r /Y — A^" 

hh,cv,x ' hh,cv,y hh,cv hh,cv,z 



A < + A< A< — 7A< 

\e—h,< ^ I hh,cv ' hh,cv,z hh,cv,x hh,cv,y 

hKcv,aa' - o I A< - 5< 

hh,cv,x ' hh,cv,y hh,cv hh,cv,z 
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A! 



~h,a 



A" + 5" 

ee,vc ' ee,vc,z 



5"' 



-I- 

X ' "ee,vc,y 



(5" 
A" 

66,1)6 



- 

66,1)6,1/ 

- 5" 



A 



/i,< 

ee,vc,aa' 



1 / A< +5< 

j 66,1)6 ' "66,1)6,2 

66,116,2) ' ee,vc,y 



6< 

"66, 116, X 

A< 

66,1)6 



"66,1)6,1/ 

• 6< 

ee,vc,z 



V. NONEQUILIBRIUM PAULI-DIRAC SPIN EQUATIONS 

The twelve transport equations for the components of the multi-band spin magnetiza- 
tion transport equations can be written is more compact form as spin-vector equations, a 
generalization of Eq. ([1]) to multi-band spin-magnetization quantum transport equations. 

The results can be summarized as four-coupled vector equations for the spin magnetization- 
distribution functions, 



_ 1 
~ 2 

1 

+ 2 



He, S^^ 



+ 



*-'e61 ^66,0 



+ i 



Bcc X S< - S< X Be 



+ 



cc cc 



+ 



■"66 ^CC,0 ^CC,0' — 'cc 



^CC ^ ^CC 



+ 



+ 



^fc^cc 



'-'cc^cc 



+ 



' — cc cc,o 



or 'z:< 

cc,o' — 'cc 



+i 



' — 'cc cc cc ' — 'cc 



1 

+ 2 



J_ 

ih 



+ 



^hh,cv'^vc 'J^-^o. 



61) 66,1)6 



hh,cv vc,o 



SI Jl 



cv,o ee,vc 



+ - 



'■'hh,cv ^ '-'ne '-'ei; ^ "66,iie 



A < ca _ or A < 

hh,CV 1)6 61) 66,116 



+ 



hh,cv vc,o cv,o ee,vc 



+i 



"hh,cv ^ '-'vc '-'cv ^ "66,i;e 
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(24) 



+i 



K? C< _ C< » 

^CC X X Byy 



1 

2 



■yr c< _ C< A « 

cc CD cc hh,cv 



+i 



"cc ^ ^CD '-'cc ^ "fe/l,CD 



1 

2 



y< qa 
cc cv 



or A < 
'^cc^hh.cv 



+ 



"cc CD,0 



''-^cc.o "hh.cv 



+i 



"CC ^ *-'cD '-'cc ^ "fth,Cl) 



1 

2 



Ar c;r>T 
hh,cv'~-'vv 



q<yrT 

CD DD 



hh,cv vv,o '-''•I' nil 



^CD.O DD 



CD DD hh,CV DD 



1 

2 



A < qrT Q*" y 

hh.cv vv ^cv vv 



'S'cD ^ 



vv 



hh,cv vv,o ^cv,o^ vv 



s:< y qrT 
"hh,cv ^ "^vv 



1 

2 



1 

+ 2 



+ 



1 

+ 2 



1 

+ 2 



^vSyc ^VC-^C 



+ 



vc.o'-^cc 



+i 

ya,T q< 

^ VV vc 



^vv X — iS^ X Be 



C'>T A a 
'-^DD ce.vc 



+i 



'^a,T q< _ q>T xa 

' — ' DD DC,0 

q>T y, sa _ 'Ba,T y, q<" 



y>T qa _ coT A < 

+i 



VV ee,vc 
•aT 



qaT w r< 

^DD ce.vc 



■B>Tqa _ qaT ?< 1 
" vv vc,o vv,o ee,vc\ 

■5>T ^ qa 



A'' S< 

ee.vc cc 



+i 



q<ya 
vc cc 



+ 



"ee.DC '^cc 



Xr q< 
ee,vc cc,o 

_ C< V '5a 
'^vc "o 



A < ca _ qr y< 
ee,vc^cc '-'vc^cc 



q< "a 
'-^ vco^cc 



A< qa _ qr 

"ee,DC ^cc,o '-^dc,o"cc 



+i 



x< sy qa _ qr Y '::;< 
"ee,DC ^ '-'cc '-'dc ^ "cc 
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1 

2 

1 

2 

1 
2 

1 
2 



+i 



+ i 



+ 



' — ' vv vv,o 



vv vv 



C< '^ct 

titi,o' — ' vv 



vv 



vv 



vv 



y< Ca or y< 

vv vv ^vv vv 



'z:< qa _ cr 'z:< 

' — ' vv vv^o vv,o' — ' vv 



+i 



X S'^.,, SI'.., X 



Ar C< _ C< A a 

ee,vc'~^cv '-^vc hh,cv 



+ 



ih 



"ee,vc'-^ cv,o '~^vc,o"hh,cv 



A< Qo- _ Q'' A< 
^ee,t)c*-^cii '-^vc'-^hh,cv 



+ 



£< ca _ or £< 

"ee,vc^cv,o ^vc,o"hh,cv 



+i 



ee,vc ^cv ^dc "hh,cv 



A. Single Conduction Band Limit 



In the single conduction band limit for the electrons, Eq. flM|) reduces to 



_ 1 
~ 2 

1 

+ 2 



1 

+ 2 



1 

+ 2 



+ 



yr C'< 
cc''-^cc 



C<ya 
'-^cc cc 



+1 

+i 



+ 



^cc,o' — 'cc 



< 

cc'—'cc 



+ 



'-'cc ' — 'cc 



^fc ^ ^cc ^cc ^ ^fc 
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Upon ignoring the effect of 5'^^o,~ the above further reduces to, 
ih 



d d 



^cc 



+ 



i?, + ReS^„5, 



< 



+ Im SL, S< S< , Im SI + 



S^^, Re 5*^^ 



+ 



2 L 



Re 'El^ X 



+ 2 



51 X Re S. 



Re X S< 



+ 2{^-XImS:, + Im x S< } . 



Further upon neglecting the effects of Re S^^,— as was done in Ref.'^, we finally obtain 
ih 



d_ _d_ 



^cc 



+ 



He + Re S^^, S^^ 



Bcc + Re El X S< - S< X Bcc + ReE 



Im H"' X S^. + X Im E''^ 



H< X Im>S,", + Im x S< 



Using the following relations, 



Then, we obtain 



ims;;. 



7, 

r. 



-Br 



eh _^ 



' 2m*c 



for electrons. 



_ 1 

~ 2 

i 

+ 2 



S^< X fC-ReS 



ReSL X 5, 



(29) 
(30) 
(31) 
(32) 



+ 2(7 x^<+ S<xl} --{E<.A + A.E<}, (33) 

which agrees with the single-band spin-vector magnetization transport equation given in 
Ref.— , within the approximation used and using Eqs. ( I29|) -( l32|) except for the presence of 
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the extra term in Eq. (133|) given by 



This terms was overlooked in casting the Re S*" single-particle Hamiltonian in Ref.= into 
its 2 X 2 spin canonical form. Equation (l28ll is the exact expression for the single-band spin 
magnetization vector transport equation. 

VI. PSEUDO-SPIN CORRELATION FUNCTIONS, S<^^ 

In addition to the twelve multi-band magnetization transport equations of the Pauli- 
Dirac spin components, we still have four more equations for the 5*^^ ^ making a total of 
16 equations. By the same token as was done for the 2x2 real-spin matrices, we can 
separate the 2x2 matrix containing S"^^ ^ into a term independent of the band indices, the 
total charge in the system, plus the pseudo-spin magnetization terms. To do this, we need 
to revert to the electron picture as discussed before. Then, the additional magnetization 
transport equations has to do with pseudo-spin magnetization by virtue of our two-band 
model, i.e., the presence of two-band discrete quantum labels v and c. 

From Eqs. f l20|) . fl2T]) . fl22|) . and f l23|) . we have the following four transport equations for 
the Pauli-Dirac spin scalar S^^ ^, where the subscripts a and (3 have the range on the set 
{v,c}, 




2 I cc cc,o cc,o cc\ ' 2 ' — ^cc cc cc ' — 'cc 




' 2 I hh,cv'-^vc,o '~-^cv,o ee,vc\ ' 2^ "hh,cv '-^vc '~-^cv 



(34) 
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+ 



+ 



K . Q< — Q< . K 

'^cc '-'cv ^Cll "^vv 



cv cc hh,cv 



'Z:< . ca C"" . X< 

cc cv cc hh.cv 



2 \_HcSf.v,o ^cv,0'^v\ 

' r\ I cc cv,o cc.o hh,cv\ 



2 L 

2 [ /i/i.ct) ^vv,o ^cv,o vv\ 2 
2 I hh,cv vv,o cv,o vv\ c2 



\y< _ Qr A < 1 

\_ cc cv,o cc,o hh,cv} 



^hh,cv '-'vv 



?< . c?rT 

hh,cv,z z,vv 



cv ' — ' 



s. 



z,cv ' — ' vv,z 



(35) 



d_ d_ 

_ 1 
~ 2 

2 L 

1 r.. 









~ Svc,oHc] + 2 


E ■ ^< 

i^'vv '-'vc 


vv,o ee,vcj 


1 

" 2 


' — ' vv 


- S''^ A< 1 - 

vv,o ee,vc] 


1 

2 


■S>T . 

' — ' vv 


_ c< f"^] 1 
,o '-^vc,o ccj ' 


1 

2 


"ee,vc 


_ qr y<~\ I 

3 '-^VC,0 Ccj ' 





-'vc '-'cc 



vc 



q>T _ ra 

vv ee,vc 

vv ecMC 



vc ' — 'cc 



Qa or 

*-^2,CC "-^ZiUC ' — 'CC,Z 



(36) 



i^vv '-'vv '-^vv '-^'v 



_l_ 1 ryr n< _ C'< 1 _|_ 1 

' 2 L vv vv,o vv,o vvj ' 2 

-L 1 ry< Q"- — 9'" y< 1 + 1 

2 L vv vv,o vv,o vvj ~ 2 

-i_ 1 TA'' '?< — S'^ A " 1 -I- - 

^ 2 l ee,vc'-^cv,o '-^vc,o hh,cv} ' 2 

-L - \A< — A< 1 -L 1 

^ 2 L eejfc ct),o '-^vc,o hh,cvi ' 2 



r _ (-<< C'< _ '^a 

ft' fi; ff ' — ' vv 



qa qr 

vv vv vv 



:< 

vv 



cr . q< _ q< . X a 

"ee,vc '-'cv '-^vc "hh,cv 



"ee,tic '-^cv '-^vc "hh,cv 



(37) 



Note that Eqs. (I24l)-(l271) and Eqs. (I34 l) -( l37l) yield 16 coupled transport equations. By 
separating S"^^ ^ into the total charge in the system plus the pseudo-spin, the multi-band spin 
magnetization quantum transport equations nonlinearly incorporates the effects of pseudo- 
spin magnetization. We have already perform the sign change of Eq. (fTOj) to obtain the 
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electron picture of Eq. (137|) above. What remains to be done for Eqs. (135|) and (136|) in 
order to revert to electron picture is to replaced all the transposed quantities, F^J with 
their equivalent expressions given in the Table 1-3. 

A. Spinless Two-Level Atom and Flat-Band Limits 

We observe that in Eqs. flM]) -f p7|) each equation contains spin-indpendent terms and 
corresponding terms which involved spin-vector dot products expressing the summation 
(i.e., the process of 'integrating out') of the real-spin degree of freedom. To gain some 
understanding into these equations, we examine Eqs. f l34|) - fl37|) by first retaining only the 
spin-independent terms, i.e., for the moment ignoring the vector dot product portions. 

We have the resulting transport equations for the scalar 5*^^ ^ involving only the spin- 
independent terms of Eqs. (1MI) - (P7I) . 




hh,cv vc,o cv,o ee,vc} ' c\ \_ hh^cv vc^o cv,o ee,vc. 



(38) 




(39) 




1 - - 1 - 

-i_ _ Ta^ Q< _i_ S'>^ A'^ 1 -I- t \A< -I- S"*^ A< 1 

' c2 I ee,vc '-^cc,o ' '~^vv,o ee,vc} ' 2 I ee,vc'-^cc,o ' '-^vv,o ee,vc} 




(40) 




(41) 
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To gain insights in these equations, as representing the pseudo-spin part of the spin 
magnetization transport equations, let us reduce these equations to a two-level flat energy 
bands or a two-level atomic system. Let us take the zero of energy in the middle of the 
band or energy gap, and let He = —Hy = We will also ignore terms arising from the 
particle self-energies, S, as well as terms involving A^. We identify the interband matrix 
elements as A^, = (^;| Hj |c), and A^^l^^^ = (c| Hi \v), so that A^^'^^^^ = A^;,^^,, and by virtue 
of the locality of space-time dependence. The above four equations reduce to the following 
expressions, 

9 1- 
d I 

'^^'0^^00,0 = 2 \P'^°^cv,o + ^\ih,cv {^vv,o ~ ^cc,o)] ) (43) 

'^^t^vco = \ [-^oS<,, + A^,,,, {S<^^ - S<J] , (44) 
(9 1- 

where we made use of the relation:— 

S'>^ = - S'< (Af\\ 

in our two-level atomic or fiat-band limiting case.— 

The pseudo-spin correlation functions are derived by first expressing the 2x2 matrix in 
the band indices into a spin-canonical form as. 



cc,o cv,o I -L 

C< C< I 2 



So^ol + So - 



^ \ Sx o ~1~ o So o Sz o 



(47) 



where. 



^Sy^O (5'^c,0 ^CV,o) ' 

c _ / c< _ c< ^ 

'-'2,0 V cc,o ^vv,o) 1 

So,0 = {S^c,0 + ^VV,o) 5 (48) 
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where we drop the superscript in the pseudo-spin correlation functions, Sj^o- We note 
that So,o represent the trace of the original 4x4 spin matrix for the two bands. Thus, So,o 
represent the total charge of the system which may vary in space and time, with the caveat 
that the background positive charge have to be subtracted from the charge represented by 
the correlation density So^o to obtain the net charge. 

Similarly, we write the 'energy'-matrix in terms of the Pauli matrix as 



H + B, B^-iBy 
B^ + iBy H-B, 



(49) 



where we defined B and H as. 



Bx — {^le,vc + ^hh,cv) ) 

iBy = (Agg „g — ^hh^ci) ) 

Bz = {He - Hv) = ^o, 

H^{H, + H,)=0. (50) 

Therefore, we have the transport equations for the pseudo-spin correlation functions, 

d 1 

^^'q^^x,o = -^[—fi^oSyco'^ ^le,vc {^cc,o ~ ^vv,o)] 
+ 2 [^oS^v,o + ^hh,cv {^vv,o ~ ^^c,o)] 

d I 

^^df^^y'" ^ 2 ^le,vc {^cc,o ~ ^vv,o)] 

+ - [—huJoScv,o ~ ^hh,cv {^vv,o ~ ^^,o)] 

2 {.^X^OjZ B^Sq^x] ) 

d 1 - 

^^^'S'z,o — 2 [^hh,cv^vc,o ~ ^cv,o^le,vc\ 

^ 2 L ee,vc'-^cv,o ' '~^vc,o hh,cvi 

— 2 {^xSo,y — BySo,x) ■ 
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In vector equation form, the pseudo-spin magnetization equation is, 



|4=iBxis„, (51) 

The above equation can also be written as 

= Bx jls (52) 

where 

Jo I — fJ-B 



mc \2 

corresponds to the magnetic moment of the pseudo-spin (pseudo-Bohr magneton), and 

—J=B 

has the units of the magnetic field, B. Thus we have realized a pseudo-spin angular momen- 
tum 1 5*0, with effective magnetic field, B, determined by the interband terms in A's and 
energy gap, Eq. (l50l) . Equations (15T1) and ( l52l) are also known as the Bloch equations. Note 
that So rotates about the 2;-axis with frequency ^ = u^r^ in the counterclockwise sense. 

Indeed, the pseudo-spin angular momentum has eigenvalues ±|/i by virtue of the two 
discrete quantum-energy labels, since from Eq. (l48ll we have, 

Sz^o = 1 for Syy^o = for the excited state, 

Sz^o = —1 for Scc,o = for the unexcited or ground state. 



B. Spin-Independent Contribution 

We now consider the full spin- independent terms of Eqs. (!38|) - (14T|) . The pseudo-spin 
correlation functions are given by Eq. f HHj) We expressed in spin-canonical form the following 
matrix in band indices. 



ImS^^ ImA^; \ I C + G Cx-iCy 



y 



(53) 



(54) 
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V< A < 

cc hh,cv 

A< y< 

ee,vc vv 



c< + c< 



y 



(55) 



Note the use of caligraphic B in Eq. ( H9|) in contrast to the use of the Greek /3 in Eq. ( l53l) . 
The calculation is tedious but the results can be expressed in vector form as 



+ 4 
1 

+ 4 



4 



1 

+ 4 



/5) 'S'o o 



I 

+ 4 



4 



^ X So - X /3 
C X + X C 



+ 4L 



C<,Re5; 



+ 



C< X ReS„"-Re5„" x 
C<,ImS^;| + |c<,Im5:„ 



+ -|C<xImS; + ImS:xC< 
where we use the following defined relation, 



(56) 



• or or or 

y,o vc,o cv^o 



(Re + i Im S^^) 



1. The Equation for the Scalar So 



We have the equation for So,o representing the total charge is given by, 

d 



ih 



1 

+ 4 

1 

+ - 



C<-^'^P Re St 



(57) 



where the superscript '^^^ means taking the dot product with respect to the pseudo-spin 
vector components. One observes that the first four terms in the the right hand side of Eq. 
(!57|) exactly correspond to the terms in the equation of the particle or charge correlation 
function in system without the spin degree of freedom. 



25 



C. Spin-Dependent Contributions 



The spin-dependent terms ignored in Sec. IVI Al represent the direct couphng of pseudo- 
spin to real (Pauh-Dirac) spins. We extract from Eqs. fl34|) - fl37|) the terms involving the 
dot products of real-spin vectors for their contributions to the pseudo-spin magnetization 
transport equations. Since the real-spin degree of freedom is essentially integrated out, what 
is left are the discrete band indices. First, we expressed all 2 x 2 quantities in the band 
indices into their pseudo-spin canonical matrix form. We make use of the expressions already 
given by Eqs. P7|) . pHj) . fj49|) . and f l50|) for the pseudo-spin canonical form of the scalar 
quantities. The vector quantities must also be expressed into their pseudo-spin canonical 
matrix form in the band indices producing dyadic tensors. 

We have for the real-spin vector dot product portion of the pseudo-spin equations, where 
the dot product can be written for convenience, using the Einstein summation convention 
for the Pauli-Dirac spin vector dot product. 



_l_ ^ fcr q< C< z:'' 1 _|_ ^ '^^ 1 

2 I — ^cc,i ^cc,i ^cc,i ' — 'cc,i\ ' 2 L' — ^cc,i ^cc,i ^cc,i ' — 'cc,i\ 

I 1 [A^ Q< _ Q< A« 1-1-1 \()< 9" _ Q'' A 

2 L ^'■h,cv,i '-^vc,i '-^cv,i "ee,vc,i} ' I hh,cv,i '-^vc,i '-'c" ^ 



< 

cv,i "ee,vc. 



(58) 



--\(n .-u-'" <?<_<?< (n .--^rT \i 

2 l\'-'cc,i ^cc,ij '-'cv,i '-'cv,i \'-^vv,i ^ vv,ij } 



1 



2 



s. 



>T1 



r\ \_^cc,i ^hh,cv,i ~^ ^hh,cv,i '~'vv,i} 



' r\ L' — 'cc,i '~^cv,i ' '~^cv,i ' — ' vv,ij 



qrT , qr s.< 1 
hh,cv,i '-^vv,i ' '~-^cc,i ^hh,cv,i\ ' 



< 



(59) 
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VC,0 



~ 2 ii^'"'"''^ ^ vv,i) ^vc,i ^vc,i {^cc,i '''' ^cc,i)] 

~ 2 L eejUCji cc,i ~ vv,i ee,vc,i\ vv,i vc,i ~ z,vc,i "cc,jj 

+ Q [^^,?;c,i 'S'cc,i + ^vv,i ^fe,vc,i\ ' (60) 



1 1 

2 I — ' ff,i vv,i ' — ' vv,ii ' 2 L'^ '>}'v,i vv,i vv,i ^ vv,i\ 

-I- 1 rA*" Q< — Q< A " 1 -1-1 rA< Q« _ A< 1 ffil ^ 

~ 2 L ee,iic,« vc,i hh,cv,ii ~ I ee,vc,i cv,i vc,i hh,cv,ii " V / 

We have the following pseudo-spin canonical matrix form for the real-spin vectors labeled 
by the discrete band indices, 



Re Sl^^, B,,,i + Re E;,^, ) 2 5^ . + . 5. _ . 

^^^le,vc,i ^^Kv,i J \Cx,i + iCy,i Ci - Cz,i 

y< A < \ T I r< A. /■< _ jc< 

cc,i hh,cv,i \ j Si ' Sz,i Sx,i Sy,i 

A< y< I 2 1 A< I ,/^< ?< _ A< 

ee,vc,i vv,i J \ Sx,i ' ^y,i ^z,i 



(62) 



(63) 



(64) 




1 I Si + Sz,i Sx,i — iSy^i 
^ \ Sx.i -l- iSyi Si — Sz 



(65) 



■'x,i I '"~'y,i '-'I >~'z,i 

We will also make use of the definition 



iSy^i = i (Re + i Im 5^^^) 

= SlcA ~ Sly , (66) 
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The results for the pseudo-spin magnetization equation from the contributions of real- 
spin vector dot product terms in Eqs. fl58|) - fl6Tl) can also be written as a pseudo-spin vector 
equation similar to the scalar-contribution given by Eq. (^E^ as, 



_ 1 

~ 4 

i 

+ 4 

1 

+ 4 



•Si, Si 



1 

+ 4 



13 i, Si 



Ci) "^M "I" 1 Ci) S, 



I 

+ 4 

1 

~ 4 



Ci ^ Si -\- Si X (^i 



1 

+ 4 



+ 



4 L 

z 

4 
1 



C<xRe5[-Re^[xC< 
C<,Im^[| + |c;<,Im5, 



+ -K^xIm5[ + Im^[xC< 



(67) 



1. The Equation for the Scalar 5, 



Similarly the contribution to the equation for So,o goes as. 



l[B.,S.]+l 



BT, s. 



i[C<,Re5[]-^z{C<,Im^[} 



1 

+ 4 



<-psp 



Re 51 



-z{c;<-^^MmS^[}, 



(68) 



where the superscripts ''^^^ indicates the process of taking the pseudo-spin vector dot products 
or the process of 'integrating out' the pseudo-spin degree of freedoms, aside from the dot 
product of the Pauli-Dirac spin vector indicated by the Einstein summation convention for 
the repeated 'i' index. 
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VII. NONEQUILIBRIUM PSEUDO-SPIN EQUATIONS 



The nonequilibrium pseudo-spin vector transport equations result by adding Eqs. fl56p 
and (Ej 



H, So 



_ 1 

" 4 
1 

+ 2 

i 



/3, Soo 



1 

+ 4 



C<,ReS^; 



+ ^{c,5^o}-^{cMms^: 



2 

+ 4 



+ ^ {C, ^0,0} -^{CMm 
C < xReSl- Re x C 



/3x So-SoX ^ 
- ^ |C X S^, + ^„ X C} + i |C< X Im^; + ImSt' x C 

1 r- -,1,- 



4 
1 

+ 4 

1 

+ 4 

i 

+ 4 



I3i, Si + — (^j^ , Re S", 
xSi-S,xE, 



2 

+ 4 



C<xReS^[-ReS^[xC;< 
- ^ {C* X + 5. X Oj + ^ {C< X ImSl + lmSlx Q 



(69) 



The corresponding absolute scalar— equation is obtained by adding Eqs. ( 1571) and ([68 

d 



ih 




(- 










\ 




1 


+ 




4 






+ 


It 




1 


+ 




4 



+ 4 



"C<-^^^ Re 5^] + ^ [C'^, So] - \ {C<-^^^ Im Si] 



(70) 



Equations (Elj) - (1271) and ( l69l) - ( 1701) constitute the nonequilibrium multi-band magne- 
tization quantum transport equations. These sixteen equations for the components of the 
spin vectors are the main results of this paper. 
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VIII. MULTI-BAND SPIN QUANTUM TRANSPORT EQUATIONS 

We summarize the main results off this paper by recasting the pertinent equations into 
more meaningful expressions. We start by rewriting the RHS of Pauh-Dirac SMQTEs into 
three groups as, 



+: 



+ 



He, S[ 
+ 



+ 



1 

+ 2 



1 

+ 2 



(^cc) ~ (^cc) 

C< (Aa \ 

cv \ ee,vcj 



< 

vc 



+ 



A < ca 
hh.cv vc 



5"" A< 

cv ee,vc 



+1 



+i 



X S'y^ S^. X I . 



+i 



cc cc 
'^hh,cv ^ '-'vc '-'cv ^ "ee 



+ i 
+ i 



s: 



X 



"hh,cv ^ '-'vc '-'cv ^ "ee,vc 



+ 



+ 



~^cc 



cc.o 



cc.o 



+ 



+ 



\ hh,cv J '-^vc,o 



cv,o \ ee,vc J 



^cc (*^cc,o) ("^cCjo) "ccj 



+ 



^< CO 

hh,cv '-^vc,o 



or s:< 
cv,o ee,vc 



(71) 



ih 



d d 



dti 



1 
2 



1 
2 



+ 



^cc'-'cv '-'cc^hh,cv 



Ar q>T _ c<yrT 
hh,cv vv cv vv 



+ 



cc'-^cv cc hh,cv 



A < orT _ or y>T 
hh,cv vv cv vv 



+i 



+i 
-i 



^cc ^ '-'cv '-'cc ^ '^hh,cv 



Bcc X X 

H< X S^, 



'-'cc ^ '^hh.cv 



hh,cv vv 



'—'1111 *— '.-111 1111 



?< w orT _ cr ■;5>T 
^ '-'vv '-'cv ^ ^ vv 



+ 



+ 



I^CC^CVM '^CVm'^VV 



+ 



'-^cc.o^^ hh.cv 



+ 



-'cc'-'cv.o 



or r< 
'-^cc.o ^hh,cv 



hh,cv vv,o cv,o" vv 



?< crT _ cr '5>T 
hh,cv vv,o cv,o^ vv 



(72) 
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1 

2 



1 

+ 2 
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1 

2< 



+ 
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vv vv *^vv vv 
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+i 

+i 



Ar Q< _ Q< /K a 

ee,vc cv vc hh,cv 



+ 



y< ca cr- y< 

vv vv *^vv vv 



A < qa _ qr A < 

ee,vc cv vc hh.cv 



+1 
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Byy X (S'.^y X Byy 
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■yu III; 



"ee.'uc ^ ^ci) '-'vc ^ "hh,cv 



+ i 
+ i 



X ca _ cr y ■:7< 



1 

+ 2 



+ 



+ 



"r C< 
" vv vv,o 



q< 'z:a 

vv^o' — ' vv 



+ 



r< ca _ cr r< 
ee,vc cv,o vc,o hh,cv 



+ 
+ 



ca 

" vv vv,o 



qr 'z:< 

vv.o' — ' vv 



sr q< _ q< X a 
ee,vc cv,o vc,o hh,cv 



(74) 



We also write the RHS of the pseudo-spin transport equation into a similar three group 
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of terms not involving the real spin, plus term invoving real spin as, 

So 



H, So 



+ 1 C,^o - 



C<,Re5„^ 
4 ,C<,Im^„^ 



+ 



13 So- So^i 13 



+ 



C< xReSl- Re Si x 



+ 



~\ {C X So + So^i C] + i {C< X Im^' + ImS^^ x C 

^ i3j X So^i So^i X i3j 

+ 1 [c<xReS^;,-ReS^;,xC;< 
1^0 X jSq j + iSo j X (^jj' 

+i{c<xims^;, + im5;, xc;< 



+ 



) ^OA 



+ 



^11 So 

■> Rs 'S'o j 

-i{C<,ImS^;,J-i{C<,Im^;^^ 



\ + 



+ 



+i /^) So.o 



C<,Re5^ 



+i {c, ^0,0} -i{c<,im5;„}_ ■ 



(75) 



Equations (ffTl) -( ffSl) and (ffOll constitute the nonequilibrium spin magnetization quantum 
transport equations of this paper, written into at least three groups of more meaningful 
terms. One note that the first group of terms of the magnetization vector equations is 
similar to that of the nonequilibrium spinless particle-correlation transport equations by 
virtue of the spin-independent transport coefficients, second group contains all the torque 
terms, and the last group generally contains coupling between spin, pseudo-spin, and charge. 
Note that that the equation for the pseudo-spin magnetization vector involves coupling to 
the total charge represented by So,o- 

IX. RETARDED GREEN'S FUNCTION AND SELF-CONSISTENCY 



Clearly, some important observables, described by the spectral correlation function, 
ImC, and occupation number correlation function, ihG^, are the governing dynamical 
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variables in quantum transport physics. Thus, strictly speaking one must also solve for the 
transport equation of the nonequilibrium retarded Green's function and self-energy in the 
presence of spin. The nonequilibrium quantum superfield formalism yields the equation for 



where v is the single-particle Hamiltonian, and all quantum labels are absorbed in the two- 
point arguments denoted by the numeral 1,2. Then similar procedure can be followed in 
bringing the above equation into a multiband matrix form. Then the submatrices can be 
expressed in the spin-canonical form as in Eq. fllSp . and one solves for the retarded spin 
correlation functions. This process will yield another 16 nonequilibrium transport equations. 

However, usually one is focused mainly in solving the spin magnetization transport equa- 
tions, and treat Eq. (1761) as separate calculations.— In most cases, the solution to Eq. ( 1761) 
are obtained by some sort of approximation and the results plugged into the spin magneti- 
zation equations as was done in the simulation of quantum transport equations for charge 
carriers in nanoelectronic devices. This is coupled with some further simplification by ignor- 
ing the effect of Re in the spin magnetization transport equations. For these reasons, 
we will no longer spend time in rigorously including Eq. fl7^ in our spin magnetization 
transport equations. 

A. Self-Consistent Electric Field and Potential 

There is also an important ingredient in the nonlinearity of the SMQTEs introduced 
through the need for the self-consistency of the potential distribution.—"— This in turn 
affects the spin-orbit coupling terms in the single particle Hamiltonian and even affects the 
corresponding many-body aspects of spin-orbit coupling.— Thus one also need to solve the 
Poisson equation for the self-consistent potential, similar to what was done in the numerical 
simulation of resonant tunneling devices. -i^ii^I Clearly then the total charge, represented by 
So^o obeying Eq. (170]) is also coupled, after correcting for the background positive charge, 
to the potential through the Poisson equation affecting the single particle Hamiltonian and 
many-body effects in a highly nonlinear self-consistent loop to be accounted for in real 
device numerical simulations, where clever approximations are eventually needed to bring 
the problem to a manageable proportion. 



G^{1,2) as^ 




(76) 
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X. PHASE-SPACE SPIN QUANTUM TRANSPORT EQUATIONS 



We will first give formal prescriptions on how to transform all the nonlocal two space- 
time-points correlation-function quantum transport equations, given in Eqs. f l2^ - f l27|) and 
( 1^ - ( 170]) into kinetic QDF transport equations defined in local phase-space points, i.e., 
in {p,q;E,t) phase space. This is achieved through the use of Buot's discrete phase-space 
transformation for condensed matter,— or his discrete formulation of quantum mechanics 
leading to the lattice version of the Weyl transformation. Then we will examine what has 
so far been attempted in the literature towards the treatment of kinetic QDF transport in 
spintronic devices. 

We observe that all of the terms entering in Eqs. - fl27|) and fIBO]) - f lTU]) involve 

commutators and anticommutators of two scalar functions and of a scalar and a vector 
functions, as well as sum and difference of cross products of vector functions. It is more 
convenient to cast everything in terms of commutators and anticommutators before making 
the necessary transformation to kinetic equations in phase space. Thus, we need to change 
the expressions involving cross product of vectors to commutator or anti-commutator, as 
the case maybe. We make use of the following identities for two vectors A and B, 

Ax B - B X A = iiiijk {Aj,Bk} , 
AxB + BxA = ke,jk [^j, B^] , 

where /j is the unit dyadic symmetric tensor or idemfactor, and iijk is the anti-symmetric 
unit tensor. As used before, the square bracket stands for the commutator and the curly 
bracket stands for anti-commutator of the two vector components separated by a comma. 
We also have terms that goes like, A x B — C x D occurring in Eqs. ( l25l) and ( l26l) . for 
example the term x S^,^ — S^^ x 5^^^^^ in Eqs. For the typical cross products 

involving four spin vectors in Eq. f l77|) . we will simply expand this as 

Ax B - C X D = UtijkAjBk - UijkCjDk (77) 
A. Transformation to {p, q, E, t) Phase-Space 

The QDF kinetic transport equations in {p, q, E, t)-spa.ce are obtained by applying the 
"lattice" Weyl transformation of the correlation function equations by using the following 
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set of identities^ (although continuum approximation is interchangeably used, this is not 
essential and we adapt the word "lattice" when referring to solid-state problems). For 
convenience, we give these here the LWtransformation of the following two space-time- 
point correlation functions: 

— + — ]f (12) ^ ih-F^ ip, q, E, t) , (78) 

where the Bloch states has the representation, {q\ p) = (j^) ^ exp | {p- q — Ei] [ Nfi? =^ 
in the continuum limit], representing a traveling wave in lattice space, with group velocity 
in a band A given by '^^^ffl . The above identities readily follow from the definition of the 
LWtransform, which in the continuum approximation can simply be written as 



d 

The second set is the LWtransform of a product of two-point correlation functions 
p = {p, —E) and q = {q, t)] in terms of "Poisson bracket operator", 

h f d^°-'> d^''^ d^"-^ d^^^ 



AB {p, q) = exp 
or in terms of integral operator 
AB {p, q) = 



2i \ dp dq dq dp 



a{p,q) b{p,q), 



(79) 



; dpdqKj{p,q;p,q) b{p,q) 
[Zirnj J 

j dpdqa{p,q)K^{p,q-p\q) , 



(80) 



where the factor ^^^^^ accounts for the proper normalization of the integration (counting of 
states in terms of unit action) in (p, q, E, t)-space, and the integral kernels are defined by 



KY{p,q]p',q') = j dudvexpi^^[{p-p')-v + {q-q')-u]j y (^p±^,qT^Y (81) 

For numerical purposes using discrete lattice points, the following expressions of Kp {p, q; p', q') 
is more preferable 

Ky {p, q] p', q') = J du dv exp | [(p — p') ' + ^ " j" y{p^u,q^v). (82) 
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Thus, we may write the LWtransform of a commutator [A, B] and an anticommutator 
{A, B} in terms of Poisson bracket differential operator, A, as 



[A, B] (p, q) = cos A [a {p, q) b {p, q)-h {p, q) a {p, q)] 

- i sin A {a {p, q) b {p, q) + b {p, q) a {p, q)} , (83) 

{A, B} {p, q) = cos A {a {p, q) b {p, q) + b {p, q) a {p, q)} 

- i sin A [a {p, q) b (p, q) - b {p, q) a {p, q)] , (84) 

where A = | • ^ — • . In terms of integral operators, we have, 

[A, B] {p, q) = J dp'dq'K^ (p, q; p, q) b {p, q) - b {p, q) (p, q; p, q) , (85) 

{A, B} (p, q) = / dp'dq'K+ (p, q; p' , q') b (p', q') + b (p', q') (p, q- p', q') . (86) 

(zvrn) J 

The above expressions simplify considerably when the LWtransforms are scalar functions. 
For this case, we have the 'lattice' Weyl transform of a commutator and anti-commutator 
of two operators, [A, B] and {A, B} , respectively, given by the following expressions 

[A, B] (p, q) = -i-g / dp'd(^K% (p, q- p', g') b (p', g') , (87) 

(ZTTft) J 

{A, B} (p, q) = j dp'dq'K^^ (p, g;p', q') b {p', q') , (88) 

where, 

(Pj Q', p', q') = J du dv exp [{p - p') ■ v + {q - q') ■ u] | 

X [ y {p + u,q - v) - y {p - u,q + v)] . (89) 

and 



(p, q;p', q') ^ J du duexp |^ [(p -p')-v + {q- q') ■ u] 

x[y{p + u,q-v)+y{p-u,q + v)]. (90) 
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XI. CONCLUDING REMARKS 



What we have accomphshed in this paper is the demonstration of the formal mathematical 
structure of nonequilibrium multiband spin-correlation transport equations. It has been 
demonstrated that spin-dependent self-energies due to many-body effects give rise to torques 
in the system. For example, the many-body effects in spin-orbit coupling has lead to the 
separation of self-energy into spin-independent part and the corresponding spin vector.— i^i^ 
Thus, it is expected that the results here will serve as a fundamental basis for constructing 
realistic transport equations, embodying various approximation schemes, for engineering 
expediency. 

The range of their validity of all simpler and manageable approximate equations can be 
assessed in the light of the present mathematical structure given in this paper This must be 
accurate enough for treating the space and time dependent spin relaxation and dephasing 
scattering mechanisms between conduction electrons, between valence holes, between elec- 
trons and holes, their coupling to pseudo-spin, and to the total charges. The total electron 
charges represented by So^o minus the positive background charge will be fed to the Poisson 
equation for self-consistency in the potential. It is expected that the accompanying physical 
interpretations will acquire deeper insights when applied to highly nonequilibrium situa- 
tions for spintronic device-performance applications, based on the fundamental structure of 
SMQTEs presented in this paper, . 

In particular, we have shown that the pseudo-spin is highly coupled to Pauli-Dirac spin 
transport equations and to the particle charge density. It is not clear how this pseudo-spin 
aspect of spin transport is described by conventional methods of Pauli-Dirac spin scattering 
physics so far employed in the literature on the theory of spin relaxation and dephasing 
mechanisms in solids. 

Acknowledgement 1 The author is grateful to Prof. R. E. S. Otadoy, R. A. Loberternos, 
and D. L. Villarin of the Theoretical and Computational Sciences and Engineering Group, 
Department of Physics, University of San Carlos, for their interest in spin quantum transport 
problems, which helps motivate the present extension of their joint previous work. 
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Appendix A: Matrix Equations for Bloch Electrons with Spin 



In the presence of spin degree of freedom, Eqs. (fT^ . (fT^ . and (fn|) become matrix 
equations. For the conduction-electron band, Eq. (E]), we have the matrix equation with 
spin indices given by. 
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where the repeated spin subscript, cr, indicates the use of Einstein summation convention 
over the spin degrees of freedom. 

For the hole band, we have the corresponding matrix equation determined from Eq. ( fT2|) 
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Appendix B: Spin-Dependent Pairing Green's Functions 



For the 'pairing' between electron and holes, we have from Eq. ( IT3|) 
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and for its nonequilibrium (reverse) 'conjugate' process,^ the matrix equation derived from 
Eq. (JUD, is. 
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